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ABSTRACT 

In this paper, we develop a novel reduced-rank space-time adaptive 
processing (STAP) algorithm based on adaptive basis function ap- 
proximation (ABFA) for airborne radar applications. The proposed 
algorithm employs the well-known framework of the side-lobe can- 
celler (SLC) structure and consists of selected sets of basis functions 
that perform dimensionality reduction and an adaptive reduced-rank 
filter. Compared to traditional reduced-rank techniques, the pro- 
posed scheme works on an instantaneous basis, selecting the best 
suited set of basis functions at each instant to minimize the squared 
error. Furthermore, we derive stochastic gradient (SG) and recursive 
least squares (RLS) algorithm for efficiently implementing the pro- 
posed ABFA scheme. Simulations for a clutter-plus-jamming sup- 
pression application show that the proposed STAP algorithm outper- 
forms the state-of-the-art reduced-rank schemes in convergence and 
tracking at significantly lower complexity. 

Index Terms — Space-time adaptive processing, Reduced-rank 
methods, Airborne radar applications. 

1. INTRODUCTION 

Airborne surveillance radar systems operate in a severe and dynamic 
interference environments, which may be composed of clutter and 
jamming. Space-time adaptive processing (STAP) employing a 
joint-domain optimization of the spatial and temporal degrees- 
of-freedom (DOFs) has become a key enabling technology for 
advanced airborne radar applications [1]. Although a significant 
increase in output signal-to-interference-plus-noise-ratio (SINR) 
compared with traditional approaches can potentially be achieved, 
the large computational complexity of full-rank STAP, which has 
order 0{M A ) due to a covariance matrix inversion operation where 
M is the dimension of the full rank filter, prohibits its application 
in real-time. Another complicated issue with full rank STAP is that 
the filter with large number of elements requires a large number of 
samples to reach its steady-state behaviour while providing slow 
convergence speed and poor tracking performance 0. A cost- 
effective technique in short-data record scenarios and, in particular, 
for systems containing a large number of elements is reduced-rank 
signal processing. 

In the comprehensive report by Ward [3| and the book by 
Klemm |4], several ad hoc reduced-rank methods were developed. 
The first statistical reduced-rank method was based on a principal- 
components (PC) decomposition of the target-free covariance matrix 
(5l . Another class of eigen-decomposition methods was based on 
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the cross-spectral metric (CSM)[6|. However, the PC and CSM 
algorithms have the problem of heavy computational load due to 
the eigen-decomposition. The family of Krylov subspace methods 
has been investigated thoroughly in the recent years. This class 
of reduced-rank adaptive filtering algorithms, including the multi- 
stage Wiener filter (MSWF) [7. 8] and auxiliary-vector filters (AVF) 
t9l llO|[Tn . projects the observation data onto a lower-dimensional 
Krylov subspace. Despite the improved convergence and tracking 
performance achieved with these methods, they are very complex to 
implement and suffer from numerical problems. The requirement 
for accurate estimation of the target-free covariance matrix is also 
an important issue which restricts the statistical reduced-rank STAP 
algorithms. 

In the present paper, we develop a reduced-rank approach to the 
STAP design utilizing an approach based on adaptive basis func- 
tion approximation (ABFA). The novel scheme employs the well- 
known framework of side-lobe canceller (SLC) [ 12 1 structure and 
consists of selected sets of basis functions that perform dimension- 
ality reduction and an adaptive reduced-rank filter. Compared to tra- 
ditional reduced-rank techniques, the proposed scheme works on an 
instantaneous basis, selecting the best suited set of basis functions at 
each instant to minimize the squared error. Furthermore, we derive 
stochastic gradient (SG) and recursive least squares (RLS) algorithm 
for efficiently implementing the proposed ABFA scheme. The re- 
sults for radar clutter and jamming suppression show a performance 
significantly better and lower complexity than existing reduced-rank 
schemes. 

This paper is organized as follows. Section[2]introduces the sig- 
nal model and in Section 3, the proposed reduced-rank STAP scheme 
is detailed, and the proposed adaptive algorithms and their complex- 
ity analysis are also presented in Section 3. Examples of the perfor- 
mance of the proposed reduced-rank STAP are provided in Section 
4 and conclusions are given in Section 5. 

2. SIGNAL MODEL 

The system under consideration is a pulsed Doppler radar residing 
on an airborne platform. The radar antenna is a uniformly spaced 
linear array antenna consisting of N elements. Radar returns are 
collected in a coherent processing interval (CPI), which is referred 
to as the 3-D radar data-cube shown in Fig. ??(a), where L denotes 
the number of samples collected to cover the range interval. The 
data is then processed at one range of interest, which corresponds to 
a slice of the CPI data-cube. This slice is a J x JV matrix which 
consists of JV x 1 spatial snapshots for J pulses at the range of in- 
terest. It is convenient to stack the matrix column-wise to form the 
M x 1, M — JN vector r(i), termed a space-time snapshot ffl|3l. 



Given a space-time snapshot, radar detection is a binary hypothesis 
problem, where hypothesis Ho corresponds to target absence and 
hypothesis Hi corresponds to target presence. The radar space- time 
snapshot is then expressed for each of the two hypotheses in the fol- 
lowing form, 



H : r = v 

Hi : r = as + v, 



(1) 



3. PROPOSED ALGORITHM 

In this section, a novel adaptive reduced-rank scheme based on adap- 
tive basis function approximation (ABFA) is proposed. Then, we 
derive the SG and RLS algorithms for efficiently implementing the 
proposed scheme. Finally, we compare the computational complex- 
ity in terms of multiplications of the proposed scheme and the exist- 
ing algorithms. 



where a is a complex gain and the vector s, which is the Mxl nor- 
malized space-time steering vector in the space-time look-direction. 
The vector v encompasses any undesired interference or noise com- 
ponent of the data including clutter c, jamming j and thermal noise 
n. These three components are assumed to be mutually uncorre- 
lated. Thus, the M x M covariance matrix R„ of the undesired 
clutter-plus-jammer-plus-noise component can be modelled as 



R„ = E{vv H } = R c + R, + R„ 



where H represents Hermitian transpose, R c = E{cc H }, Hj — 
E{jj H } and R„ = E{nn H } denote clutter, jamming and noise co- 
variance matrix respectively. In practice, the interference-plus-noise 
covariance matrix R„ is typically unknown. The common approach 
is to estimate it from the secondary data set which does not con- 
tain the signal of interest (r = v). In this context, we can refer the 
interference-plus-noise covariance matrix R„ as R. The process- 
ing chain of a general STAP radar is shown in Fig. ??(b). To detect 
the presence of targets, each range bin is processed by an adaptive 
2D beamformer (to achieve maximum output SINR) followed by a 
hypothesis test to determine the target presence or absence. The op- 
timum full-rank STAP 1 1 ] is given by 



7R _1 s 



(3) 



where 7 is a nonzero complex number. In practice, since R is un- 
known, the processor substitutes an estimation of R for R to arrive 
at the adaptive weight w. It is most common to compute the co- 
variance matrix estimate as R=-^ Ylf=i r (i) rH (i)- This approach 
is known as sample matrix inversion (SMI). Assuming Gaussian- 
distributed interference, an optimum detection statistic follows from 
the likelihood ratio test and appears as 1131 



\d(i)\ = \u, H r(i)\ f r, 



(4) 



where r\ is the detection threshold. The performance of (|4j is given 
by 



-P 2 

Pfa = exp( — ) 



Pd — I /xexp 

h 



V + p 2 ) 



(5) 



where Pfa is the probability of false alarm, Pd is the probability of 
detection, /3 = 7//\/aJ^~RaJ is a normalized detection threshold, 
Io(-) is the modified zero-order Bessel function of the first kind, 
and p equals the square-root of the peak output SINR. According 
to l[5}, Pd depends on both output SINR and the value of Pfa in 
a white noise detection scenario. If the value of Pfa is specified, 
Pd is a monotonic function of p. Thus, maximizing SINR likewise 
maximizes Pd for a fixed value of Pfa- 



3.1. Algorithm Description 

The novel scheme employs the well-known framework of SLC struc- 
ture and consists of selected sets of basis functions that perform 
dimensionality reduction and an adaptive reduced-rank filter. As 
shown in Fig. ??, we begin by using a desired target signal s (steering 
vector) as a first-stage matched filter to form a main-beam response 
do(i) = s H r(i). In our design, there are B sets of basis functions 



(2) {T 6 G C* 



, B} pre-stored as candidate projection ma- 



trices, where D is the length of reduced-rank filter. Based on Fig. ??, 
we have a set of values {yb(i)} at the output of SLC which can be 
written as a function of and the instantaneous projection matrix 
T b (i) as 



: u H (i)r(i) 



- H / - 
W (l 



TffOr' 



(6) 



where denotes the reduced-rank filter weight vector, Tb(i) de- 
notes the 6th set of basis functions and r'(i) = Br(i) is the received 
data projected in the null space of s, B = null(s) or I — is 
the block matrix. At the ith time instant, we select one matrix from 
these B sets of basis functions as the instantaneous projection matrix 
by using the decision rule 



b op t = arg min \eb(i)\ 

&£{!,... ,B} 



(V) 



where e&(i) = do(i) — yb(i)- The pre-stored matrices are con- 
structed by basis functions in the following way 



(8) 



where (f>b,d £ C denotes the dth basis function of the 6th set, 
d = 1, D and b — 1,...,B which is composed of a single "one" 
and (M ~ 1) zeros, according to the following 



b,d = [0 



1 



(9) 



where Zb, d is the number of zeros before the only element equal to 
one. We set the value of Zb,d in a deterministic way which can be 
expressed as 



Zb,d 



M 

u 



x (d- 1) + (6-1). 



(10) 



It should be remarked that other designs have been investigated and 
this structure has been adopted due to an excellent trade-off between 
performance and complexity. 

3.2. Adaptation Implementations 

Here, we describe the SG and RLS algorithms that adjust the pa- 
rameters of the reduced-rank filter based on the MSE and LS cost 
function. 



3.2.1. The SG Algorithm 

We present a low-complexity SG adaptive reduced-rank algorithm 
for efficiently implementing the proposed method. The reduced-rank 
MSE cost function is given by 

jMS E (w(i))=E{|do(i)-w K (i)T fl (i)r'(i)| 2 }. (11) 

By computing the instantaneous gradient terms of 1 11 11 with respect 
to u>(i) and introducing the positive step size n, the resulting ABFA- 
SG algorithm is thus given by 



Table 1. Computational complexity of algorithms. 



Vb{i) 


= u> H (i)Tf(i)r'(i), 


(12) 


b op t 


= arg min d (i) - t/b(i) 2 , 
6e{i,...,s} 


(13) 


r(<) 


= Tf opt r'(i), 


(14) 


e (i) 


= do(i) - ui H (i)f(i), 


(15) 


u>(i + l) 


= + fj,f(i)e* (i). 


(16) 



Algorithm 


Multiplications 


Full-Rank-SG 


2M + 1 


Full-Rank-RLS 


4M 2 + 4M 


MSWF-SG 


(£> + 1)M 2 + D 




+(3Z> + 2)M 


MSWF-RT S 


(£> + 1)M 2 + 3DM 




+2M + 4(D 2 +D) 


AVF 


D(4M 2 +4M + 1) 


+4M + 2 


ABFA-SG 


(5 + 3)13 + 2 




M 2 + M 


ABFA-RLS 


4D 2 + (B + 3)D + B 




M 2 + M 



3. 2. 2. The RLS Algorithm 

Here, we develop the RLS-version of the proposed ABFA scheme. 
To this end, let us first consider the LS design of the proposed 
scheme as 

i 

JLs(«(0) =X)^" i l*C7')-« J, WT ff (i)r'C3)| a (17) 
j=i 

where A is the forgetting factor. By computing the gradient of ( 117) 
with respect to us(i) equal to zero, we obtain 



= R x p, 



(18) 



where R = £} =1 A^f (j)r(j) H and p = £* =1 ^ 3 d* (j)r (j) 
denote the time averaged correlation matrix and cross-correlation 
vector respectively. Thus the optimal weight can be calculated re- 
cursively using the RLS algorithm, which is summarized as follows 



P(0) 
K(i) 

P(<) 
w(i) 



= J" 1 !, 



P(i - l)r(i) 



A + f ff (i)P(i- l)r(i) ' 

A _1 P(i - 1) + A _1 K(i)f H (i)P(i - 1), 

u>(z-l)+K(iK(i) 



(19) 
(20) 

(21) 
(22) 



where 5 is a small constant and I is the identity matrix. 
3.3. Complexity Analysis 

Here, we detail the computational complexity in terms of multiplica- 
tions of the proposed schemes with SG and RLS and other existing 
algorithms, namely the full-rank with SG and RLS, the MSWF with 
SG and RLS and the AVF, as shown in TableQ] where D is the length 
of the reduced-rank filter and B is the number of sets of pre-stored 
basis functions in the proposed scheme. The proposed scheme with 
both SG and RLS algorithms is much simpler than the full-rank with 
RLS, the MSWF and the AVF and slightly more complex than the 
Full-rank with SG (for D <C M). To compare the complexity of our 
proposed algorithms with other reduced-rank algorithms, let us take 
a simple example, say M — 64, D — 4 and B = 16, the number of 
multiplications for MSWF with SG is 21380, for MSWF with RLS 
is 21456 and for AVF is 66822 respectively, while our proposed al- 
gorithm has significant lower number of multiplications 4233 for SG 
and 43 16 for RLS respectively. 



Table 2. Radar System Parameters 



Parameter 



Value 



Carrier frequency (/ c ) 


450 MHz 


Transmit pattern 


Uniform 


PRF (/,.) 


300 Hz 


Platform velocity (v) 


50m/s 


Platform height (h) 


9000 m 


Clutter-to-Noise ratio (CNR) 


40 dB 


Jammer-to-Noise ratio (JNR) 


30 dB 


Elements of sensors (N) 


8 


Number of Pulses (M) 


8 



4. SIMULATIONS 



In this section we assess the proposed ABFA STAP algorithm in 
an airborne radar application. The parameters of the radar plat- 
form are shown in the Table [2] For all simulations, we assume 
the presence of a mixture of two broadband jammers at —45° and 
60° with Jammer-to-Noise-ratio (JNR) equal to 30 dB. The Clutter- 
to-Noise-ratio (CNR) is fixed at 40 dB. L = 1000 snapshots are 
simulated and all presented results are averages over 100 indepen- 
dent Monte-Carlo runs. We compare the proposed scheme with the 
full-rank filter, the MSWF and the AVF techniques for the design 
of linear receivers, where the reduced-rank filter with D co- 
efficients provides an estimate to determine whether the target is 
present or not. The signal-to-noise-ratio is set at dB. We eval- 
uate the SINR performance of our proposed algorithm against the 
number of snapshots L by comparing with all other schemes, which 
is shown in Fig. Q] The curves show an excellent performance for 
the proposed scheme, which converges much faster and has signif- 
icantly better performance than the full-rank filter, the MSWF and 
the AVF schemes at much lower complexity. In the second exper- 
iment, we calculate Pd based on ((5) by setting the false alarm at 
Pfa = 10 -10 . Fig.|2]presents Pd versus normalized SINR results 
for all schemes. The figure illustrates that the best detection perfor- 
mance is provided by the proposed scheme using RLS and SG adap- 
tive algorithms, followed by the AVF and the MSWF, and finally the 
full-rank filters. 



5. CONCLUSIONS 
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Fig. 1. SINR performance against snapshot with M ■■ 
dB, A = 0.9998, fi = 0.005. 
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Fig. 2. Probability of detection performance vs normalized SINR 
with M = 64, A = 0.9998, fi = 0.005, L = 1000 snapshots. 



In this paper, we developed a reduced-rank approach to the STAP 
design utilizing an approach based on adaptive basis function ap- 
proximation (ABFA). The novel scheme employed the well-known 
framework of the SLC structure and consists of selected sets of basis 
functions that perform dimensionality reduction and an adaptive 
reduced-rank filter. Compared to traditional reduced-rank tech- 
niques, the proposed scheme worked on an instantaneous basis, 
selecting the best suited set of basis functions at each instant to min- 
imize the squared error. Furthermore, we derived the SG and RLS 
algorithm for efficiently implementing the proposed ABFA scheme. 
The results for radar clutter and jamming suppression showed a 
performance significantly better and lower complexity than existing 
reduced-rank schemes. 
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